Andreev bound states in iron pnictide superconductors 
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Recently, Andreev bound states in iron pnictide have been proposed as an experimental probe 
to detect the relative minus sign in the s±-wave pairing. While previous theoretical investigations 
demonstrated the feasibility of the approach, the local density of states in the midgap regime is small, 
making the detection hard in experiments. We revisit this important problem from the Bogoliubov- 
de Gennes Hamiltonian on the square lattice with appropriate boundary conditions. Significant 
spectral weights in the midgap regime are spotted, leading to easy detection of the Andreev bound 
states in experiments. Peaks in the momentum-resolved local density of states appear and lead to 
enhanced quasiparticle interferences at specific momenta. We analyze the locations of these magic 
spots and propose they can be verified in experiments by the Fourier-transformed scanning tunneling 
spectroscopy. 

PACS numbers: 74.45.+C, 74.78.-w, 74.20.-z 



Discovery of superconductivity in LaFeAsOi-^FzJlJ 
and related iron pnictides [H, [U Q recently has 
reignited intense investigations on the unconventional 
superconductivity 0, As a first step toward the 

understanding of the underlying pairing mechanism, it 
is important to know the pairing symmetry in iron 
pnictides. Recent experiments on the nuclear mag- 
netic resonance^, HI, the penetration depth @, [T3, 11 



12 1 and the specific heat[13| reveal the signature of 
the nodes in the gap. However, some other experi- 
mental measurements Tj, [15 , H E3] show the oppo- 
site. In particular, the ang ular-resolved photoemission 
spectroscopy [3 0, HU HH clearly shows the presence 
of the gap at all points on the Fermi surface. Thus, the 
pairing symmetry in iron pnictides remains controversial 
at the point of writing. 

Band structure calculations (2^, [23[ show that the low- 
energy excitations mainly come from the two dimensional 
FeAs layers with Fermi surface composed of two dis- 
connected hole and electron pockets. Nesting between 
the hole and electron pockets causes the magnetic or- 
dering before the superconductivity sets in, hinting that 
the spin-fluctuation induced pairing 14,21, [24[ is plau- 
sible in these materials. Furthermore, renormalization 
group approaches 25, 26|, 27J show that the leading con- 
tender in various pairing symmetries is the extended s- 
wave (or s±-wave) pairing state where the Fermi sur- 
face is fully gapped but the gap functions for the elec- 
tron and hole pockets show opposite signs. The uncon- 
ventional s ign flip in the gap functions initiates intense 
studies[28|, [23, [30, [3l|, [32] to reveal its role in the super- 
conducting phase and also potential connection to the 
magnetic instabilities. 

Andreev bound states (ABS)|32J serves as an indi- 
rect probe for the relative minus sign in the s±-wave 
paring. In contrast to the d-wave pairing symmetry in 
cuprates[33j], ABS in iron pnictides appears at finite en- 
ergy with dispersion depending on the microscopic de- 
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FIG. 1: (Color online) A semi-infinite square lattice for iron 
pnictide with a flat edge. The surface potential Vo locates 
along the boundary sites colored in orange. The Fermi surface 
of the two-band model with tight-binding parameters ii = 
— 1, £2 = 1-3, ts = tn — 0.85 and /1 = 1.55 is shown. The 
pairing potential on the lattice and in the momentum space is 
shown here and the shaded area in the Brillouin zone denotes 
the sign change in the gap function. The blue dashed line 
stands for the transverse momentum where the hole pocket 
near the F point vanishes. 



tails. A significant drawback as an experimental probe 
is the lack for significant peak in the local density of 
states (32| in the midgap regime. In this Letter, we 
revisit this problem but treating the lattice effects in 
iron pnictides carefully with appropriate boundary con- 
ditions. Note that the work function near the sample 
boundary gives rise to an effective potential and shall be 
included. In contrast to previous results, we find signif- 
icant peaks in the momentum-resolved local density of 
states and large spectral weights in the midgap regime. 
In addition, we also predict the existence of the magic 
spots due to quasiparticle interference that can be veri- 
fied in experiments by the Fourier-transformed scanning 
tunneling spectroscopy (FT-STS)[H, 

To elaborate on the details, we start with the 
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Bogoliubov-de Gennes (BdG) Hamiltonian for the two- 
orbital model with the s±-wave pairing on the square 
lattice, 

Hbclg = X! [ t ab( r ' r ') ~ t 1 5 a,b5 r y] cl a (r)c ba (r') 

a, b r,v';<r 

+ [A; 6 (r,r')c oT (r)c 6i (r') + h.c.] , (1) 

a, 6 r,r' 

where cr labels the spins and a, b — X, Y denote the dxz 
and dyz orbitals respectively. The non-vanishing hop- 
ping amplitudes t a t on the square lattice include the in- 
traorbital hopping txx{x ± 1, y) — tyyix, y ± 1) = — ti, 
txx(x, y ± 1) = tyy (x ±l,y) = -t 2 , txx/Yv{x ± 1, y ± 
1) = ixx/Yi'^ T 1)2/ i 1) = — £3 and the interorbital 
hopping (xy(i ± l,y ± 1) = -txy(x Tl,y±l) = *4- 
The lattice constant is set to unity a = 1 through out 
the paper. With appropriate parameters, it was shown 
that the two-band model capture the essential features 

1 . The Fermi surface 



of the realistic band structure [36 
contains two electron pockets at the X and Y points and 
two hole pockets at the T and M points in the Brillouin 
zone as shown in Fig. [TJ The gap function with the s±- 
wave symmetry is described by the intraorbital pairing 
&xx/yy (x ± 1, y ± 1) = Axx/yy {x T 1, y ± 1) = A on 
the square lattice shown in the left panel of the Fig. [T] 

The translational invariance along the y direction 
makes the corresponding momentum k y a good quan- 
tum number. The open boundary ruin the translation 
invariance in the x direction but the solution can be con- 
structed by applying the generalized Bloch state [38l. 1391] - 
<& z (x) — Qz x , where z is the complex eigenvalue of the 
displacement operator, D$ z (x) = z$ z (x). The BdG 
Hamiltonian can be decomposed into different k y mo- 
mentum sectors and the eigenvalue problem is greatly 
simplified to (H z — el)$ = 0, or in explicit matrix form 



0,(2) 



where T\\ = tu/z - (^1 + ^) + hit, T22 = £22/2 - 
(M2 + A 4 ) + *222, T12 = -it u /z + it 12 z and A n = A 22 = 
Ai/z + AiZ. The effective hopping amplitudes and the 
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gap functions are in = — 1\ — 2i 3 cos k y , [i\ = -z^tuo^, 
<22 = — ^2 — 2i3COsfc y , [12 — —2tiCOsk y , t\ 2 = 2£4sin/c y , 
and Ai = 2AoCOsfc w . For given k y and e, one can solve 
for eight solutions of z with corresponding eigenvectors 
<f>. Unit modulus \z\ = 1 solutions are the usual plane 
waves and \z\ < 1 corresponds to the evanescent modes 
due to the presence of the open boundary. Note that, if 
z is a solution, we can show that 1/z is also a solution. 
Thus, we expect only four \z\ < 1 evanescent modes. 

The general solution for the evanescent modes is 
^(x) = 2^7=1 a y&y(x) an d should satisfies the bound- 
ary conditions. Due to the work function near the sam- 
ple edge, we introduce an additional Vq to incorporate 
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FIG. 2: (Color online) Dispersion of ABS with A = 0.004 
(corresponding to the realistic values of ~ 15meV for the 
superconducting gap) is shown in the upper-left panel. The 
surface potential Vb/|ti| = 0.16 is chosen here. The green 
shaded areas denote the continuum from the bulk states. The 
surface density of states is plotted in the upper-right panel 
and the momentum-solved probability density is shown in the 
lower panel. The blue dashed line stands for the transverse 
momentum where the hole pocket vanishes as in Fig. [T] 



the rise of the potential near the open boundary. After 
some algebra, boundary conditions can be casted into the 
matrix form, 



J5 c *(-l) + y*(0) =0, 



with 



B c 







Ai 


° ^ 




^22 





Ai 


Ai 





-in 


i*i2 


V 


Ai 




-til) 



V=V 



(3) 
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where 1 is a 2 x 2 identity matrix. Substituting 
the wave function 'i'(x) into Eq. |3]), the bound- 
ary conditions reduce to the homogeneous equation, 
(B C N' + VN) A = 0, where A = (01,02,03,04), N' = 
^ u ±^ 2 ,±^ 3 ,j-^ and AT = ($i,$2,$ 3 ,$ 4 ). If 

this homogeneous equation hosts nonzero solutions, the 
determinant must equal zero, det [B C N' + VN}=0. 

To obtain the full dispersion of the ABS, we need to 
solve for the generalized Bloch state for each momentum 
k y . Since the dispersion is symmetric around k y = 0, 
we only show the portion between k y = and k y = tt 
in Fig. [2] The green shaded area represents the quasi- 
particle states in the bulk. The continuum is separated 
into two sectors: the left part (with < k y < 7r/4) cor- 
responds to the Fermi surface near the T and X points 
in the Brillouin zone while the right part for the regime 
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close to the Y and M points. In the midgap regime be- 
low the continuum, the dispersive ABS is plotted as the 
red curve. Since we do not linearize the energy spectrum 
near the Fermi surface, our approach naturally eliminate 
any artifacts in the vicinity where the ABS merges into 
the continuum. This is clearly demonstrated that the 
dispersion glues back to the bottom of the continuum 
without any singularity near k y = 0, n. The most impor- 
tant feature of the dispersion is the dips located at the 
transverse momenta q\ ~ ±0.197r (near the T/X points) 
and qi — ±0.727r (near the Y/M points). These dip 
structures give rise to large spectral weights inside the 
midgap regime and also produce robust magic spots in 
quasiparticle interferences. 

First of all, the local density of states near the bound- 
ary, i.e. surface density of states (SDOS), can be com- 
puted numerically D S (E) = ±- Y, k 1*0 = 0, k y )\ 2 S[E- 
e(ky)]. The results exhibits several peaks as shown in the 
right panel of Fig. [2] However, in realistic situations, 
impurity scattering is unavoidable and those peaks liv- 
ing inside the continuum (green shaded area) hybridize 
with the plane-wave states. Thus, we only focus on 
the midgap regime where the evanescent modes are pro- 
tected. Clearly, there is significant spectral weight trans- 
fer into the midgap regime in sharp contrast to the previ- 
ous studies. Two significant peaks, arisen from the dips 
in the dispersion, can be detected in tunneling junctions 
or related experimental setups. Compared with previous 
approaches, two ingredients seem to explain the major 
differences: the band structure is kept without lineariza- 
tion (beyond the Andreev equations) and the jump of 
the work function near the open boundary is included. 
With the large spectral weight transfer and also the pro- 
nounce peaks inside the midgap regime, the ABS can be 
detected experimentally and serves as an indirect probe 
for the pairing symmetry. 

The momentum resolved probability density 
\ty(x,k v )\ 2 near the T/X points can also be com- 
puted numerically as shown in the bottom part of 
Fig. For clarity, we only show the momentum range 
< k y < 7r/4 near the T/X points. In the vicinity of 
k y = 7r/4, there is no ABS and the probability density 
vanishes. A general trend is manifest: the ABS is more 
localized near the boundary when the dispersion dives 
away from the continuum into the midgap regime. Let 
us concentrate on the peaks near the open boundary 
(red areas near x — 0) indicated by the blue and pink 
arrows first. When comparing with the dispersion, it is 
clear that the most dominant feature in the probability 
density profile occurs at the dip of the dispersion 
(indicated by the blue arrow). In fact, by integrating 
over all momenta, it is easy to check that the peaks in 
SDOS come from the dips in the dispersion. Another 
peak (indicated by the pink arrow) in the probability 
density profile arises due to the Van Hove singularity. 
Note that the blue dashed line corresponds to the 
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FIG. 3: (Color online) The dispersions of ABS with different 
surface potentials Vo/|ii] = 0.08, 0.16, 0.40, are showed as 
solid, dotted, dashed lines respectively. On the right-hand 
side, the surface density of states with Vb/|£i| = 0.08 and 0.40 
are plotted as light solid and dark dashed lines for comparison. 

transverse momentum where the hole packet near the T 
point disappears. Meanwhile, the electron pocket near 
the X point is almost empty as well. Thus, the density 
of states to form ABS is greatly enhanced and leads to 
the peak structure in the probability density. However, 
since the energy lies inside the overall continuum, it is 
unlikely to remain sharp. 

In addition to the peak features, it is also interesting to 
explore the quantum interference patterns in momentum 
resolved probability density. Close to the blue dashed 
line, stripe oscillations appear and extend into the bulk 
regime. The enhanced interference pattern comes from 
the fact that the solutions for are almost real and thus 
make the oscillation length longer and visible. On the 
other hand, near k v = regime, the solutions for Zi are 
generically complex. Linear combinations of the four so- 
lutions cancel out the quantum interferences and make 
the probability density profile more or less smooth. 

In FT-STS, one can scan the midgap regime and looks 
for magic spots in the Fourier-transformed spectroscopy 
data. Gradually increasing the energy, we will hit the 
first peak in SDOS as shown in Fig. [2] The dominant 
quasiparticle scattering comes from qi and leads to an 
enhanced interference peak at the transfer momentum 
Qi = ±2<7! ~ ±0.387r. As we further increase the energy, 
the second peak in the LDOS reigns. The quasiparti- 
cle scattering should lead to another magic spot at the 
transfer momentum Qi = ±2qi — ±1.44-7r. However, 
since there are also ABS at other momenta, we expect 
the magic spots at Qi is likely to be weaker than those 
at Q\. 

Since the ABS occurs at finite energy, it depends on 
microscopic details. Therefore, it is important to explore 
how the above predictions change when the boundary 
potential Vq varies. We show the ABS dispersions with 
several reasonable choices for the boundary potential Vb 
in Fig. [3l Clearly, the detail shape of the dispersion is 
sensitive to the boundary potential. As Vq increases, the 
dispersion dives further into the midgap regime. How- 
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ever, several generic features remain robust : the spectral 
weight transfer inside the midgap regime is significant 
and there exists a sharp peak in SDOS. Therefore, the 
predictions in previous paragraphs remain qualitatively 
correct but may suffer some quantitative corrections. 

In summary, we develop a lattice-model approach for 
Andreev bound states in iron pnictides and find it can 
be detected in experiments as an indirect probe for the 
pairing symmetry. So far, we only consider the simple 
intraorbital pairing. In recent theoretical investigations, 
interorbital and interband pairing in iron pnictides are 
widely studied and deliver many interesting and exotic 
phenomena 40l. 41 1. In principle, the methods we devel- 
oped here can be applied to the interorbital or interband 
pairing but the detail analysis remains open and requires 
further investigations. 
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Council in Taiwan through grant NSC-97-2112-M-007- 
022-MY3 and partial financial supports from the Na- 
tional Center for Theoretical Sciences in Taiwan. 
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